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ABSTRACT 

Comparisons with experimental heat transfer and 
surface pressures were made for seven turbine vane and 
blade geometries using a quasi-three-dimensional thin- 
layer Navier-Stokes analysis. Comparisons are made for 
cases with both separated and unseparated flow over 
a range of Reynolds numbers and freestream turbu- 
lence intensities. The analysis used a modified Baldwin- 
Lomax turbulent eddy viscosity model. Modifications 
were made to account for the effects of: 1) freestream 
turbulence on both transition and leading edge heat 
transfer; 2) strong favorable pressure gradients on re- 
laminarization; and 3) variable turbulent Prandtl num- 
ber on heat transfer. In addition, the effect on heat 
transfer of the near-wall model of Deissler is compared 
with the Van Driest model. 


Nomenclature 

A - Area 

Damping coefficient 
C - Chord 
Cf - Friction factor 
C p - Specific heat at constant pressure 
k - Thermal conductivity 
i - Length scale 
m - Meridional direction 
P - Pressure 

Pressure gradient, [yU^ j p)[C j j2)^^ 2 dP / ds 
Pr - Prandtl number 
Q - Dynamic head, pU 2 j 2 

Re - Unit Reynolds number, subscript gives length 
r - Radius 
s - Surface distance 

s + - Normalized surface distance, (sU € / is)y/Cj /2 

St - Stanton no., k\dT / dy\w /[{pU)uiC p (r/ N -r w )] 

T - Temperature 

Tu- Turbulence intensity 

U - Inviscid velocity 

v - Velocity in flow solution 

x - Axial direction 

Y - Pressure coefficient, Y — (P/ N — P)/Qw 
y - Perpendicular distance from surface 

- Normalized surface distance, (yU e /u)\JCj / 2 


- Intermittency factor 

S - Full boundary layer thickness 

€m - Eddy viscosity 

r\ - Direction outward from blade 

0 - Momentum thickness 

k - Von Karman’s constant, k = 0.4 
A - Momentum thickness pressure gradient 
parameter 

/x - Dynamic viscosity 
i / - Kinematic viscosity 

£ - Direction around blade 

p - Density 

<f> - Ratio of axial velocity to mean blade speed 
w - Vorticity, dv/dy 
Subscripts 

c - Edge of boundary layer 
ex - Exit 

j - End of transition 

1 - Start of transition 

in - Inlet 

m - Meridional 
s - Distance along blade 
t - Turbulent 
x - Axial 
w - Wall 

e - Tangential direction or 

- Momentum thickness 

Superscript 
7 - Total conditions 


INTRODUCTION 

Accurate predictions of turbine blade heat trans- 
fer are important, especially as turbine inlet tempera- 
ture increases. One approach to predicting blade heat 
transfer is to use a boundary layer analysis such as the 
STANS analysis of Crawford and Kays^ 1 ^ . However, this 
approach cannot be used to predict heat transfer beyond 
the point of separation. Also, because of the necessity 
of specifying a starting condition with boundary layer 
analyses, the heat transfer predictions in the leading 
edge region may be in error. An alternate approach to 
predicting turbine blade heat transfer, especially useful 
when separation occurs, is the use of Navier-Stokes anal- 
yses. Both two- and three-dimensional Navier-Stokes 



analyses have been developed* 2-6 ). The emphasis in 
these analyses has been on obtaining accurate aerody- 
namic results. This paper discusses the results of tur- 
bine blade heat transfer predictions using the quasi- 
three-dimensional, thin-layer Navier-Stokes analysis de- 
veloped by Chima* 2 ). The modifications described by 
Chima, Turkel, and Schaffer* 7 ) were used. 

The analysis developed by Chima incorporated the 
Baldwin and Lomax* 8 ) turbulence model. This model 
has several attractive features from a computational 
standpoint. However, the model does not account for 
several factors known to influence turbine blade heat 
transfer. Among these factors are: (1) the effect of 
freestream turbulence on leading edge heat transfer and 
on the location of transition, (2) the effect of bound- 
ary layer relaminarization on heat transfer, and (3) the 
effect of a variable turbulent Prandtl number. The tur- 
bulence model was modified to include models for these 
various effects. Also, the near-wall damping model of 
Deissler* 9 ) was compared with Van Driest’s model. 

The solution procedure for Navier-Stokes analyses 
differ from boundary layer analyses in several respects. 
Before alternate models can be evaluated, it is neces- 
sary to verify that other factors do not cloud the com- 
parisons. In boundary layer analyses all of the solution 
domain is confined to the boundary layer region. For 
Navier-Stokes analyses, however, only a small fraction 
of the blade-to-blade grid may be within the boundary 
layer regions. Also, boundary layer analyses are march- 
ing type analyses, while Navier-Stokes analyses require 
the solution of the entire flow field in an iterative man- 
ner. One purpose of this paper is to show results for 
different grid and iteration parameters to demonstrate 
the values needed so that comparisons are independent 
of these parameters. 

The appropriateness of the heat transfer models are 
shown by comparisons with experimental data. Com- 
parisons were made for each of the three blade rows in 
the st age- and- a-half data of Dring, Blair, Joslyn, Power, 
and Verdon* 10-13 ), the two stators of Hylton, Mihelc, 
Turner, Nealy, and York* 14 ), as well as the stator and 
rotor data of Hippensteele, Russell, and Torres* 15-16 ). 
These sources of data result in comparisons over a wide 
range of Reynolds numbers and turbulence intensities. 
The Navier-Stokes analyses have the ability to calculate 
heat transfer when there is separation. Therefore, com- 
parisons are made for cases with attached flows and for 
cases with separated flows. With respect to the thermal 
modeling there are three objectives. The first is to show 
that the accepted models for boundary layer analyses 
are appropriate for Navier-Stokes analyses as well. The 
second is to choose a model which gives reasonably good 
agreement for the cases examined. Since improved mod- 
els for turbine blade heat transfer are continually being 
developed, it is not the purpose of this paper to demon- 
strate that a particular model gives good agreement with 


data for every case. However, to demonstrate the as- 
pects of the Navier-Stokes heat transfer predictions it is 
necessary to choose a reasonably accurate model. The 
third objective is to show the effects of various assump- 
tions on turbine blade heat transfer. Since, some of these 
assumptions show better agreement than the base-case 
model, a recommendation is made for the appropriate 
prediction methodology. 


DESCRIPTION of ANALYSIS 

Basic flow analysis. The analysis used the flow 
solver developed by Chima* 2 ). The solver is for either 
inviscid or viscous quasi- three-dimensional flow, wherein 
the variation of the stream sheet thickness is accounted 
for. The variation in stream sheet thickness is deter- 
mined from a hub-to-shroud flow analysis using a pro- 
gram such as MERIDL of Katsanis* 17 ). The viscous 
analysis of reference 2 employs the thin-layer Navier- 
Stokes assumptions, wherein the streamwise derivatives 
of the viscous terms are neglected. A four-stage Runge- 
Kutta algorithm with residual smoothing and variable 
time steps is used. The equations used are given in ref- 
erence 2. The heat transfer predictions are dependent 
on an effective viscosity, /x e ff, and an effective Prandtl 
number, P£ e ff* These quantities are given by: 


A*eff — /^laminar A* turbulent — p{y "h ^A/) (l) 


*-- 1 /»+(£)/»> (2) 

The following description focuses on the calculation of 
Mturbuientj and the turbulent Prandtl number, Pr t . 

Baldwin-Lomax turbulence model. This model is a 
two-layer algebraic eddy viscosity model. The eddy vis- 
cosity is found from: 


.. _ (Mt)inner 


if y ~ YcroBsover 
if y > ycroHsover 


(3) 


where y cr ossover is given by the minimum value of y for 
which (^t)inner ” (Mt)outer» The inner and outer values 
are: 

(^)inner = p(«y(l - e _S,+ /j4 + )) M (4) 

(Pt)oxxtet = C'outerPPwAKEi'KLEB (5) 

A + = 26, and C ou t e r = 0.0269. Pwake is given by: 
PWAKE — ™l>i(y max P max > 0.25y max Vj)jp/P max ) (6) 
The quantities y max and P max are found from: 

r{y) = y | w |(l- e -«' + M + ) (7) 
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J/max is the y value for which P(y) is a maximum. Pkleb 
is given by: 

*KLEB = [l + 5.5(0.3y/y max ) 6 ] ‘ (8) 

The term vdif is given by: 

^DIF = {\J^e "h ^mlmax — (\J v> e + U m)min (§) 

Except for wakes, the second term is zero. 

When (/it)MAX in profile/^oo > 14, transition 
occurs. /Xqo is the freestream molecular viscosity. The 
turbulent Prandtl number, Pr t = 0.9. 

Turbulent Prandtl number model. As an alterna- 
tive to the constant Pr t used in the Balwin-Lomax 
model Crawford and Kays^ give Pr t as a function of 
the distance from the wall. This model is referred to as 
the Kays and Moffat model, and expresses Pr t as: 

= r 0.5 0.2 e M Pr 

r * l-Pr too v\[Prt oo 

-0.04(^^) 2 (l.O - J (g) 

The recommended value for Pr tOQ is 0.86. Figure 
1 shows the variation of Pr t with distance from the 
waD,t/ + , for a flat plate boundary layer. Pai and 
Whitelaw^ 18 ) used a linear variation in Pr t with a value 
of 1.75 at the wall and decreasing to 0.5 at the edge of 
the boundary layer. The model of Kays and Moffat was 
used because Pr t smoothly varied into the commonly 
used freestream value. 

Eddy viscosity model. In the near wall region the 
most commonly used eddy viscosity model is the Van 
Driest model. This model is incorporated into the 
Baldwin-Lomax model, and is given by equation 4. 
Other models have been proposed for the near- wall re- 
gion, White^ 19 ). One of these is Deissler’s^, where the 
eddy viscosity is given as: 

e M = n 2 vy{ 1 - ( 11 ) 

where n = 0.124. The form of equations 4 and 11 is 
very different. Away from the wall, the eddy viscosity 
becomes proportional to t/ 2 |w| in the VanDriest model. 
In Deissler’s model the ratio of eddy to molecular vis- 
cosity is proportional to the Reynolds number based on 
distance from the wall. Figure 2 shows the similarity of 
the two models for a flat plate boundary layer. 

Near wall damping model. In the Baldwin-Lomax 
model A* = 26, while in both the Kays and Moffat^ 20 ) 
and the Cebeci-Smith ( 21 ) turbulent boundary layer 
models is a function of the freestream pressure gra- 
dient, P + . The Kays and Moffatt model gives: 

A+ = 25/(C K mP + + 1) (12) 


Ckm = 30.2 for P + < 0.0, otherwise Ckm = 20.6. The 
Cebeci-Smith model gives: 

A + = 26/\/(11.8P + + 1) (13) 

Figure 3 gives the damping coefficient as a function of 
the pressure gradient for the three models. The great- 
est differences among the models occur for favorable 
pressure gradients. When A + — ► oo, the boundary 
layer completely relamin arizes. The Kays and Mof- 

fatt model shows complete relamin arization occuring 
at lower favorable pressure gradients than the Cebeci- 
Smith model. 

Crawford and Kays^ recommend that a lag equa- 
tion be used for the pressure gradient value used to de- 
termine the damping coefficient. They recommend the 
form suggested by Jones and Launder^ 22 ). 

dP + /ds+ = -(P+_P+)/C (14) 

Here P+ is equilibrium value for the pressure gradient, 
and C has a recommended value of 4000. The stream- 
wise distance between nodes in a Navier-Stokes analysis 
is much greater than the marching distance used in a 
boundary layer analysis. This can result in large point- 
to-point variations in the pressure gradient. The lag 
equation smooths out these variations, and was always 
used when was allowed to vary. 

Freestream turbulence model. The strong favorable 
freestream pressure gradients near the leading edge re- 
sults in laminar flow, even with high freestream turbu- 
lence. The freestream turbulence increases the laminar 
heat transfer. The model of Smith and Kuethe^ 23 ) was 
incorporated into the analysis to account for this effect. 
The equation for the turbulent viscosity is: 

£\f = 0.164yTu oo Z7 oo (15) 

This equation was developed for the heat transfer anal- 
ysis of cylinders in crossflow, and here refers to 
the upstream velocity. One difficulty with applying this 
model in a computational scheme is that the eddy vis- 
cosity continually increases with the distance from the 
blade. The model of Forrest^ 24 ) can be used to overcome 
this difficulty. This model gives: 

€Af = CfITUocUoo (16) 

where the mixing length, £, is defined as 

i = min{Ky, 0,0865) (17) 

Since one of the attractive features of the Baldwin- 
Lomax turbulence model is that boundary layer param- 
eters are not calculated, the model of Forrest was imple- 
mented by calculating the mixing length, i in the same 
manner as in the Baldwin-Lomax model. This was done 
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by letting 1= ny when the Baldwin-Lomax eddy viscos- 
ity was calculated from the inner formulation. When the 
outer formulation was used, the mixing length in For- 
rest’s model was either held constant at the crossover 
value or set to aero. The coefficient C F is found from: 


C F = mm(0.1875|A|/(0.25|A| + 0.01) , 0.75) (18) 


where 

V ds 


(19) 


The calculation of due to freestream turbulence was 
modified from the original form in that it was applied 
only when the flow was laminar. Otherwise, the analysis 
would give significantly increased heat transfer due to 
freestream turbulence even for fully turbulent flow. 

The momentum thickness is not calculated directly. 
Rather, Twaites’ method, as given by White^ 19 ), is used. 
The momentum thickness at a distance si along the 
blade surface is given by: 


' w “SSWJT d? * (20) 

The origin in this calculation is taken as the geometric 

stagnation point. 

Transition model. In the Baldwin-Lomax model 
transition occurs when the pressure gradient does not 
effect transition, except insofar as it affects the vortic- 
ity. Perhaps more importantly, the effect of freestream 
turbulence is also not accounted for. The choice of tran- 
sition models to be included in the calculation procedure 
was limited to those which explicitly accounted for the 
effects of both freestream turbulence and pressure gradi- 
ent. The model of Dunham^ 25 ^ was chosen. This model 
gives the start of transition Resi as : 


Re ei = (0.27 + 0.73e" 80ru ) (550 + 680/(1 - D)) (21) 


with D given as a function of the pressure gradient pa- 
rameter A: 


D = mtn(21A - lOOTu , 0.75) (22) 

Where A is given by equation 19. Comparisons for a 
range of A and Tu values showed Dunham’s model to 
give lower Ree% values than the model of Seyb^ 26 ). Dun- 
ham’s model gave nearly the same values when com- 
pared with the model of Abu-Ghannam and Shaw^ 27 ^. 
Dunham’s model was chosen for implementation over 
the other models because the pressure gradient enters in 
a linear fashion, and is less likely to result in oscillations 
in the transition location as the calculations progress. 

Even the approx imate calculation of the momen- 
tum thicknesses can cause difficulties, especially near 
the stagnation point. The edge velocity, U e) is calculated 
from the local static and freestream total pressures. If 


stagnation occurs other than at the geometric stagna- 
tion point, U € will — ► 0 and 9 will — + oo. Consequently, 
transition would occur at this point. To prevent this, 
transition was not allowed to occur within a small frac- 
tion, (about 10%), of the the surface distance from the 
stagnation point. 

Transition length is determ ned from Dhawan and 
Narashima^ 28 ) model. The end cf transition is given by: 

Rest = Re si + cResi' S (23) 

The coefficient, c, has a value of 5.0 for a transition 
length between 25% and 75% intermittency. For the full 
intermittency range between 0 and 99 %, the coefficient, 
c, increases by a factor of 3.36. The intermittency factor 
7 1 of Abu-Ghannam and Shaw^ 27 ) is used: 

— 4.6b((Rcs —i2esi)/(J?esf — -R^si t . 

7t = 1 - e V / (24) 


CALCULATION PROCEDURE 

The GRAPE code of Sorenson^ 28 ) was used to gen- 
erate C-grids for the analysis. The locations of the inter- 
section of grid lines with the inner and outer boundaries 
was controlled to give minimum shear in regions of in- 
terest. For example, it was often desired to have near 
normal grids on the rearward portion of the suction sur- 
face. However, because of the large amount of turning, 
a near-normal grid could not be generated over the en- 
tire blade surface and still maintain grid-line periodicity 
on the outer boundary. Fortunately, the presence of 
sheared grid lines did not appear to adversely affect the 
flow solutions. 

The thick trailing edges made it necessary to use 
trailing edge cusps. Without a cusp, there was a large 
loss in total pressure in the trailing edge region, even 
for the inviscid calculations. Inviscid calculations were 
used to determine the shape of the cusp. The cusp shape 
was adjusted to give no net loading over the cusp region 
beyond the physical end of the turbine blade. 

After the Euler solution was obtained, the Navier- 
Stokes solution was calculated for the same blade ge- 
ometry, but with a t/ + « 1 for the first node from the 
wall. For the Euler calculation the first node spacing 
wall increased about an order of magnitude to avoid an 
unacceptably large total pressure loss. 

Steady state was reached in 3000 to 5000 time cy- 
cles. Residuals were reduced at least two orders of mag- 
nitude. Figure 4 shows the convergence history as well 
as pressure distributions and heat transfer results for 
a typical case. Reductions in residuals by two orders 
of magnitude were needed for converged heat transfer 
results, but fewer cycles were needed for pressure. 

Figure 5 shows converged results for three differ- 
ent grids. Grids of 145X46, 284X46, and 145X31 were 
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used. The coarser grid in the £ direction, 145X46, gave 
satisfactory agreement with the finer grid, 284X46, in 
this direction. However, the results were sensitive to 
the number of grid lines in the r} direction. The number 
of grid lines in the £ direction was determined by the de- 
sire to resolve the surface pressure distribution. When 
strong surface pressure gradients were encountered, a 
large number of £ grid lines were used. 

Figure 6 shows heat transfer results as a func- 
tion of the distance of the first node from the wall, 
y£ y expressed in terms of EF . This value is calcu- 
lated prior to the viscous calculation, where y^ EF = 
(yi U e ju) y/Cf~j2. The value of C f is obtained from a fiat 
plate correlation: Cf/ 2 = 0.0287(sC/ c /z/)~° 2 , to give: 
^ref = 0.17yii2e°‘ 9 /s° 1 (25) 

The largest value of the unit Reynolds number, typically 
the one based on exit conditions, gives the largest value 
for y^EF* This, in turn, gives a conservative estimate 
of y*. yj^EF 13 on ly slightly affected by the choice of s, 
and s can be taken as the chord. When a —+ 0, near the 
stagnation point, so does U ey and yJ E f remains finite. 
Figure 6 shows that y+ values approximately equal to 
one are sufficient. Since y + is related to the normalized 
temperature by the Prandtl number, sensitivity to y + 
reflects the sensitivity to normalized temperature. 

The Navier-Stokes and Euler solutions took approx- 
imately 250 and 200 seconds on a Cray X-MP for 3000 
time cycles using a 145 by 46 grid. 

Predictions were made with a specified wall tem- 
perature approximately equal to the average of the ex- 
perimental wall temperatures. Varying the wall tem- 
perature showed no change in the heat transfer as long 
as the analytic gas-to-wall temperature difference was 
equal to or greater than the experimental one. It was 
felt that it was appropriate to compare heat transfer 
calculated assuming constant wall temperature with ex- 
perimental data generated from nominally constant heat 
flux conditions. In the experiments there was some lat- 
eral conduction within the blade surface, whi:h changes 
the boundary condition from one of constant flux to one 
with more slowly varying temperature. If calculations 
are made using specified wall flux, comparisons of heat 
transfer coefficients may not adequately verify a code’s 
ability to predict heat transfer. If the Reynolds num- 
ber is sufficiently high, and the specified heat flux suf- 
ficiently low, the gas temperatures adjacent to the wall 
will not significantly change from the no heat flux case. 
Under these conditions, comparisons of heat transfer co- 
efficients will show good agreement, independent of how 
far the first grid point is from the wall. 

The temperature gradient at the wall accounted for 
the non-orth on ality of the grid lines, and was evaluated 
using a three point formula: 

dT _ £ m r?g/r - fjmte/r (— 37\ -f 4T 2 — T 3 ) . . 

d V \A?rn + {ve/r ) 2 2 


COMPARISONS with DATA 

In the comparisons that follow, a prediction which 
gave good agreement for a number of cases is given 
by the solid curve. This model incorporates Dunham’s 
transition model, and the Smith and Kuethe freestream 
turbulence augmentation effect. In addition to this 
curve, other curves are shown in each figure to illus- 
trate the effects of various parameters on blade heat 
transfer. A further discussion of the appropriate heat 
transfer model will occur after the results for all of the 
comparisons have been presented. Conclusions based on 
one set of data may not be appropriate for another set, 
and the appropriate model is a best fit model to all of 
the data. 

Stator results 

Low speed cascade tests. Figure 7 compares pre- 
dicted and experimental results for the stator of Dring 
et al. Figure 7a shows reasonably good agreement be- 
tween the analysis and the experimental pressure distri- 
bution. The pressures are slightly higher on the rear- 
ward portion of the suction surface. The pressure gradi- 
ent, though, is accurately predicted. The pressure gra- 
dient is more important than the small differences in 
levels in predicting the heat transfer. 

Figure 7b shows heat transfer comparisons at de- 
sign Reynolds number, and low freestream turbulence. 
The transition model is very sensitive to the freestream 
turbulence at low turbulence intensities. This can be 
seen by comparing the prediction at a turbulence level 
of 1.5 % with the prediction of no freestream turbulence. 
The Baldwin-Lomax transition model gave a reasonably 
accurate prediction of suction surface transition for this 
case. It was observed that this model predicted tran- 
sition when the suction surface pressure gradient be- 
came positive, and was somewhat accurate only when 
the freestream turbulence level was low. 

Figure 7c shows the effect of a variation in the 
freestream turbulence on the transition location for the 
same Reynolds number as in figure 7b, but at higher 
freestream turbulence. One curve assumes Tu = Tu OQ) 
and the other assumes Tu = TuooUoo/Ue. Since the suc- 
tion surface pressure is near a minimum at transition, 
the second assumption results in transition occuring fur- 
ther back on the suction surface. Dunham sugested that 
the turbulence intensity be taken as the average of two 
values, which are closely approximated by these two val- 
ues. Consequently, these two curves reprent the range 
of turbulence intensities for use in the transition model. 
Also shown in this figure is a fully turbulent prediction. 

Stanton numbers at an intermediate Reynolds num- 
ber for two Tuqo values are shown in parts d and e. Fig- 
ure 7d shows that the shorter transition length agrees 
better with the rapid rise in heat transfer. Part e shows 
the effect of a variable Pr t to be relatively small, even 
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though figure 1 shows Pr* to be twice as large at the 
wall for the variable Pr t assumption. Part f shows that 
the effects of both Deissler’s near wall turbulence model 
and the Cebeci-Smith relaminarization model are also 
relatively small for the lowest Reynolds number with an 
upstream turbulence grid. 

Transonic cascade tests. Figure 8 gives comparisons 
for the Mark II stator of Hylton et al. Figure 8a shows 
that because of the rapid rise in suction surface pressure 
near midchord, a grid with 271 circumferential lines was 
needed to accurately predict surface pressures. Due to 
the large adverse suction surface pressure gradient the 
flow separated, and boundary layer analyses could not 
be used downstream of the separation. Figure 8b shows 
heat transfer comparisons for the two grid densities and 
for a case where the measured wall temperatures were 
used for the boundary conditions. This case, which used 
the 271 x 53 grid, shows that the variation in wall tem- 
perature did not account for the overshoot in predicted 
heat transfer. Even though the overshoot was large, it 
was very narrow. It appears that the thermal modeling 
is incorrect immediately after separation, or perhaps, 
there was some smoothing of the experimental data due 
to conduction. The 146 x 53 grid gives results in better 
agreement with the heat transfer data, though not with 
the surface pressure data. 

Figure 8c shows that Baldwin-Lomax transition cri- 
terion is less satisfactory than Dunham’s model for pres- 
sure surface transition using the 146 x 53 grid. Figure 
8d shows that neglecting the effect of freestream tur- 
bulence at Tu = 8.3% underpredicts the heat transfer 
in the leading edge region. Figure 8e shows that us- 
ing a variable Pr t increased the predicted heat transfer 
slightly, and did not agree better with the data for ei- 
ther pressure or suction surfaces. Figures 8b through 8e 
show reasonably good agreement between the predicted 
and experimental suction surface heat transfer down- 
stream of separation. This shows the usefulness of the 
Navier-Stokes approach for turbine blade heat transfer. 

Comparisons are shown in figure 9 for similar condi- 
tions as shown in figure 8 for the Mark II stator. Figure 
9a shows good agreement in surface pressures for the 
C3X stator of Hylton et al. There is no indication of 
separation for the higher solidity C3X stator. Figure 9b 
shows an accurate heat transfer prediction for the case 
with low Re and pressure ratio of 0.582, except on the 
rear portion of the pressure surface. This figure also 
shows that the effects of Deissler’s near wall turbulence 
model is relatively small. A prediction with Tu^ re- 
duced to 5% is also shown, and is in somewhat better 
agreement with the data. This is shown becausethe tur- 
bulence may have decayed somewhat between the mea- 
suring station and the leading edge of the blade. Figure 
9c gives comparisons at high Re, where freestream tur- 
bulence effects on heat transfer are expected to be large. 
Forrest’s model for freestream turbulence effects is very 


close to the data in the leading edge region. This model 
agrees with the data better than the Smith and Kuethe 
model. Calculations with Forrest’s model showed only 
small changes in heat transfer when the turbulent vis- 
cosity was calculated with a mixing length limited to the 
crossover value, or even neglected in the outer region. 

Figure 9d shows good agreement with the heat 
transfer data at the low Re and pressure ratio of 0.490. 
The figure also shows that the A+ model had little effect 
on the heat transfer. Figure 9e shows a fully turbulent 
calculation to agree better with the data than a model 
with transition for the high Re and low pressure ratio 
case. The improved agreement occurs in the leading 
edge region, and just prior to suction surface transition. 

Linear cascade tests. Figure 10a compares surface 
pressures for the vane data of Hippensteele et al. Figure 
10b compares heat transfer results at design Reynolds 
number and low turbulence intensity for two transition 
lengths and for the Baldwin-Lomax transition criteria. 
In contrast to figure 7d, the longer transition length 
agrees better with the data. However, either transi- 
tion length in Dunham’s model agrees better than the 
Baldwin-Lomax criterion. Figure 10c shows that the 
freestream turbulence model of Forrest agrees with the 
heat transfer data better than the model of Smith and 
Kuethe at hig h Tu. Figur es lOd shows a longer transi- 
tion length than the data for the low Re- low Tu case. 
Since even the suction surface transitioned near the rear 
of the vane, any turbulence generated within the passage 
would move transition forward. Increasing Tu to 3% in 
the prediction gave the correct suction surface transi- 
tion. Figure lOe shows that, even with high Tu, suction 
surface transition is delayed to the near midchord as a 
result of the low Reynolds number. 

Rotor results 

Large scale rotor. Figure 11 compares pressure dis- 
tributions and heat transfer results for three different 
values of flow coefficient, <j>. Changing 4> changes the 
inlet flow angle. There is good agreement in pressure 
distributions between the analysis and data for all three 
<f> values. Changing inlet flow angle gives different pres- 
sure distributions, with different boundary layer growth 
rates, which in turn affect the blade heat transfer. When 
<f> = 0.68, the flow is more tangential, resulting in a 
greater adverse pressure gradient along the pressure sur- 
face. When <f> = 0.96, the flow angle is reduced from the 
design case. This in turn gives a less favorable gradient 
along the suction surface. 

Heat transfer comparisons are shown for three <f> 
values at design Reynolds number, and for <f> = 0.78 at 
a low Re, Since the test occurred in a rotating rig, the 
rotor, being behind the stator, was always in a high tur- 
bulence environment. The effect of augmenting the vis- 
cosity to account for freestream turbulence is shown in 
figure lib. Figure lid shows that a constant value of A + 
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gives the best prediction on the suction surface. It was 
observed that for high Tu cases neglecting relaminar- 
ization often gave the best agreement with the data for 
the pressure surface. Since high freestream turbulence 
causes laminar boundary layers to become turbulent, it 
may also inhibit relaminarization of turbulent boundary 
layers. Figure Ilf shows that a fully turbulent bound- 
ary layer calculation gives a reasonably accurate predic- 
tion for the entire blade surface. However, because this 
model does not account for the level of freestream tur- 
bulence, the leading edge region heat transfer may in 
some cases be in error. The heat transfer prediction for 
the low Reynolds number case in figure llg is accurate, 
except for the prediction of transition. The underpre- 
diction of pressure surface heat transfer again suggests 
that relaminarization is not appropriate in a high Tu 
environment. 

Large scale exit vane. The exit vane of the tests of 
Dring et al. had a shape similar to the rotor. Therefore, 
comparisons for this blade row are included with other 
rotor results. Figure 12 compares surface pressure and 
heat transfer predictions. The analysis underpredicts 
the heat transfer even for the turbulent assumption on 
the rearward portion of the suction surface. While only 
one prediction is shown in figure 12, several other cases 
were run with other assumptions for Pr ti A + , and the 
near wall damping model. No case gave heat transfer 
results that agreed with the data on the rear portion of 
the suction surface. It appears that the thermal mod- 
eling is in error, or perhaps there are three-dimensional 
flow effects. The fully turbulent boundary layer predic- 
tion shown by Dring et al. also underpredicted this heat 
transfer. 

Linear cascade rotor. The rotor tested by Hippen- 
steele, Russell, and Torres showed pressure surface sep- 
aration. Figure 13a gives comparisons of surface pres- 
sures. While there is good reasonably good agreement 
between the analysis and experimental surface pres- 
sures, the analysis shows a more stronger suction surface 
adverse pressure gradient than might be inferred from 
the data. 

The tests were done in a low freestream turbu- 
lence environment. Consequently, transition occurred 
towards the rear of the blade on the suction surface. 
Figure 13b shows an early prediction of suction surface 
transition at the highest Reynolds number. This was 
probably the result of an overly strong suction surface 
adverse pressure gradient. At the lower Reynolds num- 
bers, figures 13b and c, Rt$ was not large enough to 
initiate transition at this point on the blade. This figure 
shows only a small effect on suction surface transition 
when the freestream turbulence is assumed to vary with 
the local freestream velocity, 

On the pressure surface separation occurred near 
the leading edge. At all Reynolds numbers the pres- 
sure surface heat transfer reached a second peak as a 


result of separation. Because this occurred close to the 
leading edge, 0, calculated from equation 20, was too 
small to initiate transition in this region. Best agree- 
ment with the pressure surface data was achieved when 
the Navier-Stokes analysis was used to indicate transi- 
tion. Once a negative wall friction was calculated, the 
flow was assumed to be fully turbulent. This assump- 
tion gave the location of a transition like increase in heat 
transfer only slightly downstream of separation. The 
higher heat transfer along most of the pressure surface 
compared with the suction surface was predicted. The 
analysis predicted an increase in heat transfer at separa- 
tion. The increase was from a lower laminar level, and, 
in contrast to the experimental data, the heat transfer 
did not reach the leading edge level. The results in figure 
8 showed that a denser grid changed the heat transfer 
predictions at separation. For the rotor case, however, 
the heat transfer did not change with a denser grid. 

Choice of Heat Transfer Model 

The results in figures 7 through 13 s,tow that the 
choice of heat transfer model is somewhat dependent on 
the flow field environment. When both the Reynolds 
number and Tu are high, a fully turbulent model gives 
accurate predictions. When the Reynolds number is low, 
even if Tu is high, transition must be accounted for. 
Dunham’s transition model gives accurate predictions, 
and the actual level of freestream turbulence, even if it 
is low, should be used. When the turbulence level is 
low, but not known, a value of at least 2% should be 
used. When the turbulence level is low, and especially 
if the Reynolds number is also low, a negative skin fric- 
tion should be used to switch to fully turbulent flow. 
This should not be done in the immediate region of the 
stagnation point so as to avoid spurious transition. At 
stagnation the sign of the shear stress is indeterminate. 
The transition criterion in the Baldwin-Lomax model is 
less accurate, principally because it fails to account for 
freestream turbulence. The data comparisons did not 
indicate a strong preference for either the 0 to 99% in- 
termittency model or the 25% and 75% one. Both gave 
reasonably accurate results. 

The choice of relaminarization model also was 
dependent on the freestream turbulence environment. 
With high Tu it did not appear appropriate to in- 
crease when the pressure gradient was strongly fa- 
vorable. Since high freestream turbulence causes a lam- 
inar boundary to become turbulent, it also may inhibit 
relaminarization of a turbulent bounary layer. Other- 
wise the Kays and Moffat model appeared slightly fa- 
vorable, perhaps only from a computational standpoint, 
to the Cebeci-Smith model. 

Not only did Deissler’s near wall turbulence model 
give similar results to the widely used VanDriest model 
for a flat plate boundary layer, the results were nearly 
the same for the heat transfer after separation occurred. 
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Using a variation in Pr t instead of a constant value 
did not result in significantly different heat transfer re- 
sults. However, there was some indication that at low 
Reynolds numbers a variable Pr t would give higher heat 
transfer rates. 

The model of Forrest for the heat transfer augmen- 
tation in laminar flow was in better agreement with the 
data. If the model of Smith and Kuethe is used, best 
agreement with the data is achieved when the coeffi- 
cient used in their model is reduced by about 50%. A 
reduction of this magnitude is consistent with the ob- 
servations made by O’Brien and VanFossen^ 30 ^. 

CONCLUSIONS 

Comparisons with experimental data showed good 
agreement in terms of surface pressures for the several 
different turbine blades. When the surface pressure 
changed abruptly, grids of high density were required. 
Otherwise, a grid size of 145 x 46 was adequate. 

Reductions in residuals of at least two orders of 
magnitude were required in order to achieve converged 
heat transfer predictions. Convergence was slower for 
the pressure surface heat transfer than for the suction 
surface heat transfer. Use of the GRAPE code to gen- 
erate grids proved satisfactory. 

The ability of the thin-layer Navier-Stokes analysis 
to accurately predict turbine blade heat transfer down- 
stream of separations demonstrates the utility of this 
approach. After separation, a fully turbulent predic- 
tion agrees well with the data. The degree of agreement 
between the analysis and the experimental data is pri- 
marily determined by the choice of model for transition. 
The Navier-Stokes results should be used to change the 
calculation procedure from laminar to turbulent when 
separation is indicated. 
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FIGURE 1. - VARIATION OF TURBULENT PRANDTL NUMBER WITH 
DISTANCE FROM THE WALL FOR FLOW OVER A FLAT PLATE. 



FIGURE 2. - COMPARISON OF VAN DRIEST AND DEISSLER EDDY 
VISCOSITY MODEL IN NEAR WALL REGION. 



FIGURE 3. - COMPARISON OF KAYS AND MOFFAT, CEBECI-SMITH, 
AND BALDWIN-LOMAX DAMPING COEFFICIENTS MODELS AS A 
FUNCTION OF PRESSURE GRADIENT, 
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FIGURE N. - CONVERGENCE HISTORY FOR A SAMPLE CASE. 
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FIGURE 5. - COMPARISON OF STANTON NUMBER DISTRIBUTION 
FOR DIFFERENT GRIDS. 
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FIGURE 6. - EFFECT OF FIRST NODE SPACING ON SURFACE 
HEAT TRANSFER. 
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(e) STANTON NUMBER AT 62 PERCENT OF DESIGN Re, TURBU- 
LENCE GRID. 
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FIGURE 7. - PRESSURE DISTRIBUTION AND STANTON NUMBERS FOR STATOR - 1 OF DRING ET AL. 
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FIGURE 8. - PRESSURE DISTRIBUTION AND STANTON NUMBERS FOR MARK II STATOR OF HYLTON ET AL. 
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FIGURE 9. - PRESSURE DISTRIBUTION AND STANTON NUMBERS FOR C3X STATOR OF HYLTON ET AL. 
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(d) STANTON NUMBER AT 22 PERCENT OF DESIGN Re, NO 
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FIGURE 10. - PRESSURE DISTRIBUTION AND STANTON NUMBERS FOR STATOR OF HIPPENSTEELE ET AL. 
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FIGURE 11. - PRESSURE DISTRIBUTION AND STANTON NUMBERS FOR ROTOR OF DRING ET AL. 
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FIGURE 11. - CONCLUDED. 
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FIGURE 12. - PRESSURE DISTRIBUTION AND STANTON NUMBERS 
FOR SECOND STATOR OF DRING ET AL. 


18 



NASA 

National Aeronautics and 
Space Administration 


Report Documentation Page 


1. Report No. 

NASA TM- 102496 


2. Government Accession No. 


4. Title and Subtitle 

Navier-Stokes Analysis of Turbine Blade Heat Transfer 


3. Recipient’s Catalog No. 


5. Report Date 


6. Performing Organization Code 


7. Author(s) 
R.J. Boyle 


8. Performing Organization Report No. 

E-5219 


10. Work Unit No. 
505-62-21 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


11. Contract or Grant No. 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


13. Type of Report and Period Covered 
Technical Memorandum 


14. Sponsoring Agency Code 


15. Supplementary Notes 

Prepared for the 35th International Gas Turbine and Aeroengine Congress and Exposition sponsored by the 
American Society of Mechanical Engineers, Brussels, Belgium, June 11-14, 1990. 


16. Abstract 

Comparisons with experimental heat transfer and surface pressures were made for seven turbine vane and blade 
geometries using a quasi-three-dimensional thin-layer Navier-Stokes analysis. Comparisons are made for cases 
with both separated and unseparated flow over a range of Reynolds numbers and freestream turbulence inten- 
sities. The analysis used a modified Baldwin-Lomax turbulent eddy viscosity mode. Modifications were made to 
account for the effects of: 1) freestream turbulence on both transition and leading edge heat transfer; "2) strong 
favorable pressure gradients on relaminarization; and 3) variable turbulent Prandtl number on heat transfer. In 
addition, the effect on heat transfer of the near-wall model of Deissler is compared with the Van Driest model. 


17. Key Words (Suggested by Author(s)) 
Turbine 
Heat transfer 
Navier-Stokes 
Blade 


18. Distribution Statement 

Unclassified — Unlimited 
Subject Category 34 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of pages 

22. Price* 

Unclassified 

Unclassified 

20 

A03 


NASA FORM 1626 OCT 86 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 



OWomAi; 

OF POOR 


PAGE is 

Quality 



(a) SURFACE PRESSURES. 



SURFACE DISTANCE, s/s tot 
(C) STANTON NUMBER AT 55 PERCENT OF DESIGN Re. 






SURFACE DISTANCE, s/s to{ 


SURFACE DISTANCE, S/S tot 


(b) STANTON NUMBER AT DESIGN Re, NO TURBULENCE GRID. 


<d> STANTON NUMBER AT 37 PERCENT OF DESIGN Re. 


FIGURE 13. - PRESSURE DISTRIBUTION AND STANTON NUMBERS FOR ROTOR OF HIPPENSTEELE ET AL. 


19 



